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ABSTRACT 

Multi-frequency radio observations of the radio bridge of a powerful classi- 
cal double radio source can be used to determine: the beam power of the jets 
emanating from the AGN; the total time the source will actively produce jets 
that power large-scale radio emission; the thermal pressure of the medium in 
the vicinity of the radio source; and the total mass, including dark matter, of 
the galaxy or cluster of galaxies traced by the ambient gas that surrounds the 
radio source. The theoretical constructs that allow a determination of each of 
these quantities using radio observations are presented and discussed. Empirical 
determinations of each of these quantities are obtained and analyzed for 22 radio 
sources; Cygnus A is one of the sources in the sample, but it is NOT used, or 
needed, to normalize any of these quantities. 

A sample of 14 radio galaxies and 8 radio loud quasars with redshifts between 
zero and two was studied in detail; each source has enough radio information to 
be able to determine each of the physical parameters listed above. The beam 
power was determined for each beam of each AGN (so there are two numbers for 
each source); these AGN are highly symmetric in terms of beam powers. Typical 
beam powers are about 10^^ erg s~^. No strong correlation is seen between the 
beam power and the core - hot spot separation, which suggests that the beam 
power is roughly constant over the lifetime of a source. The beam power increases 
with redshift, which is significant after excluding correlations between the radio 
power and redshift. The relationship between beam power and radio power is 
not well constrained by the current data. 

The characteristic or total time the AGN will actively produce a collimated 
outflow is estimated. Typical total lifetimes are ~ (10^ to 10®) years. Total source 
lifetimes decrease with redshift. This decrease in total lifetime with increasing 
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redshift can explain the decrease in the average source size (hot spot - hot spot 
separation) with redshift. Thus, high-redshift sources are smaller because they 
have shorter lifetimes; note that higher redshift sources grow more rapidly than 
low redshift sources. 

A new method of estimating the thermal pressure of the ambient gas in the 
vicinity of a powerful classical double radio source is presented. This new esti- 
mate is independent of synchrotron and inverse Compton aging arguments, and 
depends only upon the properties of the radio lobe and the shape of the radio 
bridge. A detailed radio map of the radio bridge at a single radio frequency 
can be used to estimate the thermal pressure of the ambient gas. Thermal pres- 
sures on the order of 10"^'^ dyne cm~^, typical of gas in low- redshift clusters 
of galaxies, are found for the environments of the sources studied here. It is 
shown that appreciable amounts of cosmic microwave background diminution 
(the Sunyaev-Zel'dovich effect) are expected from many of these clusters. This 
could be detected at high frequency where the emission from the radio sources is 
weak. 

The total gravitational mass of the host cluster of galaxies is estimated using 
the composite pressure profile and the equation of hydrostatic equilibrium for 
cluster gas. Total masses, and mass-density profiles, similar to those of low- 
redshift clusters of galaxies are obtained. Thus, some clusters of galaxies, or cores 
of clusters, exist at redshifts of one to two. The redshift evolution of the cluster 
mass is not well determined at present. The current data do not indicate any 
negative evolution of the cluster mass, which may have implications for models 
of evolution of structure in the universe. 

Subject headings: galaxies: active — quasars: general — radio continuum: 
galaxies 
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1. Introduction 

Classical double radio sources, also known as FRII radio sources (Fanaroff & Riley, 
1974), have been extensively studied. Numerous surveys of radio sources as well as 
detailed studies of individual sources have provided a wealth of radio data on many 
sources. Theoretical models for the sources developed and discussed by various authors 
{e.g. Blandford & Rees, 1974; Scheuer, 1974, 1982; Begelman, Blandford, & Rees, 1984; 
Begelman & Cioffi, 1989; Eilek & Shore 1989; Gopal-Krishna & Wiita 1991; Daly 1990, 
2000) lend much understanding to the physics of the sources. It is believed that FRII 
sources are powered by highly coUimated outflows from an active galactic nucleus (AGN), 
and that the radio emission from these sources is the result of an interaction between the 
beam, or jet, and the ambient medium. As a result, careful study of FRII sources can 
yield very useful information about the FRII sources, their environments, and their central 
engines (e.g. Daly 2000). 

Daly (1994, 1995) presents a model for powerful extended FRII sources and showed 
how key parameters of an FRII source and its environment, such as the beam power and 
the ambient gas density, can be estimated using multi-frequency radio data. Results for a 
sample of powerful 3CR radio sources are also presented by the same author (Daly 1994, 
1995). A larger sample was compiled and studied in detail; results from these studies have 
been presented in a series of papers (Wellman & Daly 1996a,b; Wellman 1997; Wellman, 
Daly, & Wan 1997a, b [hereafter WDW97a,b]; Wan & Daly 1998a,b; Guerra & Daly 1996, 
1998; Guerra 1997; Daly, Guerra, & Wan 1998; Guerra, Wan, & Daly 1998; Wan 1998; 
Daly 2000). For example, detailed results on the density and temperature of the ambient 
medium of the FRII sources are presented by WDW97a,b, and the redshift evolution of 
the characteristic size of an FRII source and its application for cosmology are presented by 
Guerra & Daly (1996, 1998), and Guerra, Daly, & Wan (2000). 
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This paper presents results on several key parameters of the FRII sources and their 
environments that were not presented and discussed in previous papers. These parameters 
include the luminosity in directed kinetic energy of the jet, also known as the beam power, 
Lj, the total time the AGN will produce the coUimated outflow, i*, the thermal pressure of 
the ambient medium, Pth, the effect of the hot gas on the microwave background radiation, 
and the total gravitational mass of the host cluster (we find that the radio sources lie at the 
centers of gas rich clusters of galaxies). The way that these parameters can be estimated 
from radio data, and the empirical results for the sample, are presented here. 

The paper is structured as follows. A brief description of the sample and data for the 
sample are presented in §2. Discussions on the jet luminosity and timescale of jet activity 
are given in §3, the thermal pressure of the ambient medium and estimates of the effect 
of the hot gas on the microwave background radiation are presented in §4, and estimating 
the total gravitational mass of the the host cluster of the radio source is discussed in §5. 
Each of these sections are divided into two subsections, with the first subsection presenting 
theory and the second one presenting empirical results for the sample. Further discussions 
of the results and a summary of the paper are presented in §6. 

Values of Bubble's constant Hg = 100 h km s~^Mpc~^ and the de/acceleration 
parameter qo=0 (an open empty universe with zero cosmological constant) are used to 
estimate all the parameters in this study. Different choices of cosmological parameters, 
within reasonable limits, give results that are very similar to, and consistent with, those 
presented here. 
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2. Sample Description and Data 

The sample used in this study has been described in detail in several previous papers 
{e.g. WDW97a,b; Guerra & Daly 1998). The readers are referred to these papers for a 
full description of the sample; a brief description of it follows. All the FRII sources in 
this sample are very powerful 3CR sources, having 178 MHz radio power greater than 
3 X 10^^ h~'^ W Hz~^sr~"^; such powerful FRII sources are referred to as FRIIb sources 
(Daly 2000), or Type 1 FRII sources (Leahy & Williams 1984). They are drawn from the 
samples of Leahy, Muxlow, & Stephens (1989; hereafter LMS89), which contains sources 
with large angular sizes, and Liu, Pooley, & Riley (1992; hereafter LPR92), which contains 
sources with smaller angular extent. The final sample contains 27 radio lobes from 14 
FRIIb radio galaxies, and 14 lobes from 8 FRIIb radio loud quasars. These sources have 
redshifts ranging from zero to 2, and projected core-hot spot separations between 25h~^ 
kpc and 200/i~^ kpc. Note that the questions of projection effects and radio power selection 
effects have been addressed in great detail by Wan & Daly (1998a,b). 

WDW97a,b used the radio maps from LMS89 and LPR92 to estimate the width of the 
radio bridge behind the radio hot spot, a^, the non-thermal pressure of the radio lobe that 
drives the forward shock front, Pl, and the rate at which the bridge is lengthening, referred 
to as the lobe propagation velocity vl- The radio information may also be used to estimate 
the beam power and total time the source will be active, as described in §3. 

In order to estimate the Mach number of lobe advance, high quality radio maps 
which image large portions of the radio bridge are required; the radio bridge is defined as 
the radio emitting region that lies between the radio hot spot and the origin of the host 
galaxy or quasar. There are 16 bridges in the sample for which we have a detection of 
the Mach number of lobe advance, which includes 13 radio galaxy bridges and 3 radio 
loud quasar bridges. The Mach number may be combined with the lobe propagation 
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velocity vl to estimate the temperature of the ambient gas (WDW97a). Thus, the ambient 
gas temperature and pressure may be estimated for these sources. Lower bounds on the 
Mach number, and hence upper bounds on the ambient gas temperature and pressure, 
are available for the 14 other bridges, including 8 galaxy bridges and 6 quasar bridges, as 
described in WDW97a,b. The maps of the sources do not have sufficient dynamic range 
to allow detailed analyses of their bridge structure, and thus do not have estimates on the 
Mach number of lobe advance, nor on the ambient gas pressure and cluster mass. As a 
result, the sample of sources with estimates of the ambient gas pressure and cluster mass is 
rather small. 

Results for radio-loud quasars and radio galaxies are analyzed separately when there 
are enough data points, as is the case in the study of the beam power. As noted above, 
the number of sources with thermal pressure and cluster mass estimates is rather small (16 
lobes with detections and 14 with bounds). Thus it is not practical to separate quasars and 
galaxies in the study of these parameters. 

Cygnus A (3C 405) is the only source with a redshift ~ in our sample, and in some 
cases, results seem to depend on whether or not this source is included. Thus, subsamples 
both with and without Cygnus A are analyzed, and the results included here. Note that 
Cygnus A is not used to normalize any of the quantities studied here, with the exception of 
the total source lifetime. Note, however, that an identical normalization of the total source 
lifetime is obtained when Cygnus A is excluded from the analysis (see the companion paper 
by Guerra, Daly, & Wan 2000). Thus, Cygnus A is not needed to normalize any of the 
quantities studied here. 
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3. Luminosity in Directed Kinetic Energy of the Jet 

3.1. Theory 

The jet or beam of an FRII radio source refers to the colhmated outflow which carries 
energy from the AGN. When the jet impinges upon the external medium, a strong shock 
is formed and the kinetic energy of the jet is deposited in the vicinity of the shock front, 
which is marked by the radio hot spot. In principle, one can calculate the beam power, 
dE/dt, carried by the jet by studying the propagation of the hot spots. However, in practice 
it is difficult to use hot spot properties to estimate the beam power, since the hot spot is 
generally not resolved, and the position and properties of radio hot spots vary on short time 
scales (Laing 1989; CariUi, Perley, and Dreher 1988; Black et al. 1992.). 

The model presented and discussed by (Daly 1994, 1995) bypasses these difficulties by 
using properties of the more stable radio lobes to estimate the luminosity in directed kinetic 
energy, or the beam power, Lj. This model is a natural progression following the work of 
Blandford & Rees (1974), Scheuer (1974), Eilek & Shore (1989), Gopal-Krishna & Wiita 
(1991), Rawhngs & Saunders (1991), and many others. The method and results presented 
here are complementary to those of Eilek & Shore (1989) in the sense that Eilek & Shore 
(1989) considered the case of a non-supersonic expansion of the radio emitting region, while 
we consider the supersonic expansion of this region, and the shape of the radio emitting 
region considered here is different from that considered by Eilek & Shore. The results 
presented here are also complementary to those of Gopal-Krishna & Wiita (1991), who 
consider likely explanations for the decrease in FRII radio source size with redshift, and who 
present a consistent picture in which the increase in radio luminosity and decrease in radio 
source size with redshift can be understood. The work presented here differs from that of 
Gopal-Krishna & Wiita (1991) in that neither an efficiency nor an evolution of the ambient 
gas density with redshift are adopted. Instead, we use the radio properties of a source 
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to solve directly for the beam power of the source, as well as several other parameters, 
including the ambient gas density. The work presented here on the beam power is similar 
in spirit to that of Rawlings & Saunders (1991). The key difference is that here the radio 
bridge is assumed to be cylindrically symmetric, and thus has a volume that depends on the 
length of the bridge times the cross-sectional area of the bridge, while Rawhngs & Saunders 
(1991) seem to take a bridge volume that is proportional to the cube of the bridge length. 

Observations and numerical simulations indicate that radio lobes of powerful extended 
radio sources propagate supersonically relative to the ambient medium (Alexander & Leahy 
1987; Prestage & Peacock 1988; Cox, Gull, & Scheuer 1991; Daly 1994; WDW97a,b). Leahy 
(1990) shows that the properties of the radio lobe can be used to estimate the rate of energy 
input. In §3.4.4.3, Leahy (1990) shows that the jet power may be written (3P + P)V/t, 
since the energy density of a relativistic plasma is 3 times the pressure P; note that the 
time rate of change of the volume may be written: V/t — 7ra| vl since the forward region 
of the source grows with a velocity vl and has cross-sectional area 7ra|,. His equation may 
be rewritten as: 

Lj = At^oIvlPl- (1) 

Here is the cross-sectional radius of the radio lobe, vl is the lobe propagation velocity, 
and the FRII source is assumed to be cylindrically symmetric about the radio axis. Using 
typical units, Lj can be expressed as 

L, = 3.6 X lO^^erg s"^ { ( , (2) 

^ ^ Vkpcy V c / VlO-^°dynecm-V 

All three parameters used to estimate Lj in eq. (2) can be estimated from radio 
data (see, for example, Daly 1995, or WDW97a,b). The lobe radius can be measured 
from the radio map (see WDW97a,b for detail). The lobe propagation velocity Vl can 
be estimated using multi-frequency observations of radio lobes and bridges, by modeling 
the spectral aging of relativistic electrons due to synchrotron and inverse Compton losses 
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(c/. WDW97a,b; Wan & Daly 1998a,b; Myers & Spangler 1985; Alexander & Leahy 1987; 
LMS89; LPR92; Carilli et al. 1991; Perley & Taylor 1991). The lobe pressure Pl is given by 



where Bmin is the magnetic field strength in the radio lobe under the minimum-energy 
condition (c/. Miley 1980; Pacholczyk 1970), b is the ratio of the true magnetic field 
strength B to the minimum-energy magnetic field strength: B — hBjnim and the equation 
is in cgs units. WDW97b estimate h to be about 0.25, with a source-to-source dispersion 
less than about 15%. This value of h is consistent with values obtained by other authors 
(c/. Carilli et al. 1991; Perley & Taylor 1991; Feigelson et al. 1995), and is used throughout 
this study. 

The radio spectral index of an FRIZ source is used in the spectral aging analysis to 
estimate the lobe propagation velocity vl- It is shown in WDW97b that there is a clear 
correlation between the radio spectral indices of the sources in this sample and redshift. 
The spectral index a can be expressed as a oc (1 + zY . The best fitted value of s is 
0.8 ± 0.2 with a reduced of 1.7. Such a correlation could be caused by intrinsic curvature 
in the initial electron energy spectrum or effects of inverse Compton cooling on the hot 
spot spectral index (see WDW97b). Thus, the observed spectral index may not be the 
appropriate index to used in spectral aging analysis. This lead WDW97b to also consider 
correcting the observed spectral indices to zero redshift using the empirical correlation 
between a and (1 -|- 2;) for the sample. Results obtained with and without this correction 
are included here. 

The beam power Lj can be used to estimate the total time that the AGN will be active, 
and will produce highly collimated outfiows. Following Daly (1994, 1995) and Guerra & 
Daly (1996, 1998), the total lifetime of an outfiow, defined as i*, is related to the energy 




(3) 
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extraction rate Lj by a power law: 

U = CLJ^*/\ (4) 

where (3^ is a parameter that has imphcations for models of energy extraction from the 
central engine. As described in detail in a companion paper (see Guerra, Daly, & Wan 
2000), eq. (4) is equivalent to assuming that the source will have an average size close to the 
average size of the full population located at a similar redshift, and that the average rate 
of growth of the source will be close to the current rate of growth of the source. Thus, we 
could have written ~< D > /vl where < D > is the average size of the full population 
of FRIIb sources at similar redshift, and vz, is the rate of growth of the source under 
consideration. Since this is precisely the manner in which the value of /3* is obtained, it is 
an equivalent approach. Note that the value of is, fortuitously, independent of choice of 
cosmological parameters, and it is independent of the way the relation is normalized (see 
the companion paper of Guerra, Daly, & Wan 2000). 

The current best estimate of is 1.75 ± 0.25 (see the compansion paper of Guerra, 
Daly, & Wan 2000). (Note that the used here is not related to /So used in §4 and §5 for 
the King density profile.) The characteristic size of an FRII source is related to its lifetime 
as = viti, = CvlLj , where C is the normalization in eq. (4). This normalization 
factor is chosen so that at z ^ 0, the characteristic size of Cygnus A matches the observed 
average lobe-lobe size for sources at this redshift. That is, {vl)CL^^*^^ for one side of 
Cygnus A is added to {vl)CLJ^''^^ on the other side of Cygnus A, and this is set equal to 
2(< D >z=o) = 2 X (68 ± 14) kpc. This equation is then solved for the normalization factor 
C; the uncertainty on C is indicated in Table 1. This normalization is then used in equation 
(4) to obtain an estimate of the total time that the AGN will be producing large-scale jets. 
Note that a nearly identical normalization is obtained when Cygnus A is excluded from the 
data analysis, as shown in a compansion paper (Guerra, Daly, & Wan 2000). Note also that 
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the constraints on (3 and on cosmological parameters are independent of this normahzation, 
as described in Guerra, Daly, & Wan (2000). The value for estimated using independent 
information from each side of the source should of course be equal, and are generally very 
close in value. The value adopted for the beam power, Lj, of Cygnus A used above is 
hsted in Table 1, vl for Cygnus A is obtained from Table 1 from WDW97b, and a value of 
< D >,=o = (68 ± 14) kpc is adopted from Table 3 of Guerra, Daly & Wan (2000). 

The results on do not depend on the normalization adopted here; as shown in a 
companion paper, nearly identical results obtain when Cygnus A is excluded from the 
analysis. And, Cygnus A is not used to normalize any other quantity. Thus, none of the 
results presented in this paper rely upon a normalization based on Cygnus A. 

3.2. Empirical Results 

The beam power, Lj, has been estimated using eq. (2) for the lobe propagation velocity 
estimated using the observed radio spectral index, and the redshift-corrected radio spectral 
index; both are listed in Table 1. Typical values of Lj are 10^^ erg s ^. 

No strong correlation is found between the beam power Lj and the linear size of the 
source, represented by the core-hot spot separation r. Figure 1 is a log-log plot of Lj vs. r, 
where the observed radio spectral index a is used. The figure of Lj vs. r after applying the 
redshift-correction on a is almost identical to Figure 1, and is thus not shown here. 

For radio-loud quasars, the best-fitted line in Fig. 1 has a slope of about zero. A 
Spearman rank correlation analysis shows that the correlation coefficient between Lj 
and r is —0.03(7.72%) when no redshift-correction on the spectral index is applied, and 
—0.01(1.78%) when the spectral index correction is applied. Here the number in parentheses 
is the significance of the correlation. These results suggest that the correlation between Lj 
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and r for radio-loud quasars is insignificant, either with or without the a — z correction. 
That is, Lj is independent of r for the radio-loud quasars. 

For radio galaxies, Lj appears to increase slightly with r in Fig. 1, and results from the 
Spearman rank correlation analysis hint a marginally significant correlation between Lj and 
r. The correlation coefficient is about 0.3, with a significance level of about 90%, which 
holds with or without the a — z correction. However, note that the best-fitted line for radio 
galaxies in Fig. 1 has a slope only 2a away from zero, and that a correlation coefficient of 
about 0.3 is rather weak. Thus it appears that Lj is at most weakly dependent on r for 
radio galaxies. These results for radio-loud quasars and radio galaxies are consistent with 
the assumption that Lj is roughly constant over a source's lifetime. 

The relationship between Lj and redshift and radio power (P^) has been studied in 
detail in Wan & Daly (1998a). Increases of Lj with z and with Pj. are observed (see 
Figures 2 and 3). However, the radio power of the sources in the sample also increases with 
redshift since the sources used in this study all come from the flux-limited 3CR survey. As 
a result, which of the Lj — z and Lj — correlations is more significant, and the role of 
radio power selection effects, needs to be considered carefully. Wan & Daly (1998a) used 
two parameter fitting and partial rank analysis to investigate this. The two parameter fit 
Lj (X PiT'si^ + zY^ yields the following Uz values: Uz — 1.45 ± 0.32(0.60) for all galaxies, 
Uz — 3.83 ± 0.93(1.73) for all galaxies except for Cygnus A, and Uz — 9.56 ± 1.37(3.91) for 
all radio-loud quasars. 

The number in parenthesis is the uncertainty on times y^^, which includes the 
effect of a reduced (Xr) ^^^^ greater than one for the fit. This notation will be used 
throughout this paper. 

The values of are more that 2(7 away from zero, for both radio galaxies and 
radio-loud quasars. This means that the Lj — z correlation is unlikely to be caused purely 
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by radio power selection effects. Tfiat is, a real increase of Lj with redshift exists, though 
the magnitude of this redshift evolution is not well determined. 

The Lj — Pr correlation, on the other hand, is rather poorly constrained. Results from 
the one parameter fit suggest Lj goes roughly proportional to Pns (Figure 3). The two 
parameter fit Lj oc ^1*78(1 + 2;)""^, which takes into account the redshift dependence of radio 
power, gives values of Up with large uncertainties: Up — 0.76 ± 0.11(0.20) for all galaxies. 
Up = 0.25 ± 0.22(0.41) for all galaxies except for Cygnus A, and Up = -0.54 ± 0.29(0.82) for 
all radio-loud quasars. Thus, the exact relationship between Lj and radio power is not clear 
at the present time. Note, however, large amounts of scatter are seen in the Lj vs. Pns 
plot (see Figure 3), which is also reflected in the large reduced of the one parameter fit. 
For the same beam power, the radio power can vary by about a factor of ten. This suggests 
that radio power is not an accurate measure of the beam power, and that the apparent 
correlation between Lj and Pr may result from the fact that both increase with redshift. 

The properties of FRIIb radio sources suggest that the total lifetime of a source 
decreases systematically with redshift. This follows because the average size of FRIIb 
sources decreases with redshift, while the rate of growth of the sources increases with 
redshift. Another way to state this is: as Lj increases with redshift, the total lifetime of an 
outflow, i*, estimated using eq. (4), decreases with redshift. The values of the full lifetime 
of the source U are listed in Table 1. Typical values of range from about 10^ to 10^ 
years. Figure 4 plots as a function of redshift. It is clear that decreases with z, which 
can explain the fact that the average size of powerful extended 3CR sources decreases with 
redshift for z > 0.3 (see Guerra & Daly 1998), while the lobe propagation velocity vl 
increases with redshift. Clearly ii vl increases with redshift, and the mean source size 
decreases with redshift, then must decrease with redshift since ~ Vit^. The sources 
at high redshift produces more powerful jets (with larger Lj) for a shorter period of time 
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(with smaller i*), which results in smaller average sizes than low-redshift sources. 



4. Thermal Pressure of the Ambient Gas 
4.1. Theory 

The studies of WDW97a and b show that both the ambient gas density and temperature 
of an FRIIb source can be estimated from radio data. The thermal pressure of the ambient 
gas Pth is obviously proportional to the product of the density Ua and the temperature T. 
Interestingly, as shown below, although the lobe propagation velocity enters into Ua and T, 
it cancels out in the product, so the thermal pressure of the ambient gas can be estimated 
using single frequency radio data if the radio emission from the bridge is mapped over a 
large enough region of the bridge. This offers a completely new method to estimate the 
thermal pressure of the ambient gas surrounding high-redshift radio sources using only 
single frequency radio data. 

It can be complicated and time consuming to estimate the ambient gas density in the 
vicinity of classical double radio sources using X-ray observations, especially for sources at 
high redshift. An attractive alternative method of estimating Ua in the vicinity of a radio 
source is to use the radio data. Given that the radio lobe represents a strong shock front, 
the lobe pressure Pl, the ambient gas density pa, and the lobe propagation velocity are 
related: the strong shock jump conditions imply that Pl ~ 0.7 5 Pavj^ holds, where Pl is the 
non-thermal pressure inside the radio lobe, and Pa is the mass density of the ambient gas. 
The electron number density of the ambient gas Ua can be estimated using 

lApnipVi 

where p, is the mean molecular weight of the gas in AMU, with a value of 0.63 when solar 
abundances are assumed, and nip is the proton mass. 
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The ambient gas temperature can be estimated by careful studies of the radio bridge, 
as discussed in detail in WDW97a. A brief description of the theory is as follows. When 
the non-thermal pressure of the radio bridge is much greater than the ambient gas pressure, 
the bridge undergoes supersonic expansion. This lateral expansion can be treated as a blast 
wave and the expansion velocity is determined by ram pressure confinement. This causes 
the width of the bridge, a, to vary with time approximately as t^/^ (e.^'. Begelman & Cioffi 
1989; Daly 1990; WDW97a). Since the lobe propagation velocity vl is roughly constant 
during a given source's lifetime {e.g. LPR92; WDW97a,b), the lobe front will advance a 
distance of vit during t as the bridge expands. This means that at a distance x from the 
hot spot, the bridge expansion time is t = x/vl- Thus the shape of the radio bridge roughly 
follows a square root law, with the bridge width a{x) oc x^/^. This behavior was in fact 
found by WDW97a. As the bridge expands, the pressure inside decreases until it becomes 
comparable to the ambient gas pressure. When this occurs, the lateral expansion velocity 
becomes sonic, causing a break in the functional form of a{x). The width of the bridge 
starts to deviate from the square root law and becomes roughly constant as equilibrium is 
being reached. At the point where such a break occurs, the lateral expansion velocity is 
approximately sonic. Thus the sound speed of the ambient gas (q) can be estimated using 

■VL, (6) 

b 

where Vs\j^ denotes the lateral expansion velocity at the break. 

The Mach number of lobe propagation is defined as M = vl/cs- Given the expression 
for Cs in eq. (6), the Mach number of lobe propagation is: 



1 da 
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dx 
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^'^ dt 


b 


dx 




dx 



c, V dx 



-1 



^. (7) 



where is the width of the bridge at the break, and Xb is the position of the break relative 
to the hot spot (see WDW97a). It can be seen from eq. (7) that M depends only on the 
geometrical shape of the radio bridge. 
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The ambient gas temperature T can be expressed as: 



^^«!oc«iM-^ (8) 



where 7 is the ratio of specific heats of the ambient gas. 

Multi-frequency radio maps are required in order to use eqs. (5) and (8) to estimate 
the ambient gas density and temperature of an FRII source, since both Ua and T depend on 
vl- The product Ua T is independent of vl, so vl is not needed to estimate the 
thermal pressure of the ambient gas. The ambient gas pressure in the vicinity of the 
radio lobe is simply 

Pth = {ria + ni)kT (X (^^^ vlM-^ oc (9) 

with Pl estimated using eq. (3) and M estimated using eq. (7). Here rii is the number 
densities of ions, and for a gas with solar abundances, the electron density is Ua — 1.21ni. 
Substituting in the numerical constants in eqs. (5) and (8), the thermal pressure is given by 

Ptf^ = 0.8^, (10) 
for a gas with a specific heat ratio 7 = 5/3. 

The thermal pressure estimated using eq. (10) depends only on the lobe pressure Pl 
and the geometrically determined Mach number M, both of which can be obtained from 
single frequency observations of the radio lobe and bridge. No spectral aging analysis is 
needed in order to estimate Pth, since the dependences of Ua and T on vl cancel. 
As a result, Pth is also independent of whether or not the redshift-correction on the radio 
spectral index is applied. 

At low-redshift. X-ray data can often be used to estimate Ua and T, and hence Pth- 
However, high-resolution X-ray data are often difficult to obtain for sources at high redshift. 
This new method of estimating Pth, using single frequency radio data, offers an attractive 
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alternative. It provides a powerful tool to study the environments of powerful classical 
double radio sources, especially those at high-redshift. Since it has been shown that these 
sources are in cluster-like gaseous environments, the sources may be used to study evolution 
of gas in clusters of galaxies (as discussed by Daly 1995; WDW97a,b; Daly 2000). 



The value of Pth, the thermal pressure of the ambient gas in the vicinity of the radio 
lobe, is hsted in Table 1. Most sources in the sample have Pth on the order of (10~^^ to 
10^^°) /i^/'' dyne cm~^, which is typical of gas in low-redshift clusters of galaxies. This 
is consistent with results obtained by WDW97(a,b), who find cluster-like density and 
temperature for the ambient gas of the FRII sources in the sample. Note that the thermal 
pressures obtained here do not depend on a spectral aging analysis, whereas density and 
temperature estimates do. 

The X-cluster around Cygnus has been observed by Arnaud et. al. (1984) and Carilli, 
Perley, & Harris (1994). The thermal pressure of the ambient gas near its radio lobe is 
estimated to be about 10~^°dyne cm~^ (see CariUi, Perley, & Harris 1994). This estimate 
from X-ray data is consistent with the thermal pressure estimated here for Cygnus A (3C 
405) using the new method just described (see Table 1). 

The composite pressure profile, Pth as a function of the core-hot spot separation r, is 
shown in Figure 5. It can be seen that Pth decreases with r. A negative pressure gradient is 
expected for an isothermal gas distribution that follows the King density profile, which is 
suggested to be the ambient gas distribution for sources in our sample (WDW97a,b). For 
this gas distribution, the thermal pressure decreases with r as 



4.2. Empirical Results 




(11) 
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where Pth,c is the thermal pressure at the center of the cluster, and Tc is cluster core radius. 

The pressure gradient seen Figure 5 appears to be consistent with that expected for an 
cluster gas distribution that is isothermal and follows a King density profile, with a cluster 
core radius Tc of ~ (50 to 150) kpc. The sources in our sample have values of r ranging 
from (25 to 250) kpc. Within this radius range, the pressure profile given by eq. (11) 
has an average slope of (~— 1.4 to — 0.7), for Tc values of (50 to 150 kpc), assuming 
Po — 2/3. These expected slopes are consistent with the best-fit slopes of ~— 1.0± 0.3 in 
Figure 5. 

There is some hint of a redshift evolution of the cluster core radius from the data. 
Figure 6 plots Pth vs. r in two redshift bins, z < 0.9 and z > 0.9. The division ai z = 0.9 
is chosen so that the two redshift bins cover about the same range in redshift and contain 
about the same number of data points. It appeeirs that Pth decreases less rapidly with 
r at high redshift than at low redshift, suggesting a larger at high redshift. 
The best-fit slope of ~ —1.4 in the low-redshift bin is consistent with Tc ~ 50 kpc 
for the pressure profile given by eq. (11), whereas the best-fit slope of about —0.4 in the 
high-redshift bin is consistent with Tc ~ 250 kpc. This result is still preliminary since 
only sources with detections of Pth are included in the fits. Better estimates of Pth for 
sources currently with only upper bounds on Pth will help to better determine whether Tc is 
evolving with redshift. 

The results obtained above are consistent with results obtained by WDW97b. They 
study the ambient gas density profile and find that Tc ~ 50 kpc when high- and low- 
redshift source are considered together. They also find that that data are consistent with 
a constant core gas mass model where Tc increases with redshift roughly as Tc oc (1 -|- zy^^. 
Note that results from the study of Pth do not depend on spectral aging analysis, whereas 
those from the study of Ua do. It is thus encouraging that consistent results are obtained 
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from the two studies. 

For the following analysis, we use both a non-evolving core radius of 50 kpc, and 
an evolving core radius of 50 (1 + zY'^h'^ kpc, with a focus on the latter. 

The thermal pressure at the center of the cluster {Pth,c) can be estimated by scaling the 
pressure in the vicinity of the radio lobe (Pth) to the cluster center. This follows because 
the studies of Daly (1995) and WDW97a,b indicate that these radio sources are located 
at the centers of clusters of galaxies, so the core-hot spot separation can be used as an 
estimate of distance from the cluster center. The current data on T, ria, and Pth are all 
consistent with the gas distribution being isothermal and following a King profile. Thus 
we use the pressure profile for such a gas distribution, as given by eq. (11), to estimate the 
central thermal pressure. The central pressure is simply 



The values of Pth,c for the sources in our sample are hsted in Table 1, where a value of 
Po — 2/3 has been assumed. Most sources have values of Pth,c around 10~^° /i^/^ dyne cm~^, 
which is rather typical of gas pressure in the core regions of low-redshift galaxies clusters. 

Results on the redshift evolution of the central gas pressure Pth,c are somewhat 
uncertain because Pth,c seems to correlate with both redshift and radio power when either 
correlation is considered separately (see Figures 7 and 8) , and it is not clear at present which 
correlation is more significant. For the non-evolving core radius model, a two-parameter 
fit of Pth,c oc (1 + z)''^P^{s yields n, = 0.29 ± 0.49(0.67) and Up = 0.66 ± 0.19(0.26) when 
all sources are included, and = 1.51 ± 2.11(2.99) and Up = 0.34 ± 0.58(0.82) when 
Cygnus A is not included. For the evolving core radius model, the two-parameter fit gives 
= -1.22 ± 0.49(0.65) and Up = 0.88 ± 0.19(0.25) when all sources are included, and 
= 1.56 ±2.11(2.82) and Up = 0.15 ± 0.58(0.78) when Cygnus A is not included. Note that 
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Cygnus A appears to lie in a cooling-flow region while most of the other sources are not in 
coohng-flow regions (see WDW97a; §5). Thus fits with and without Cygnus A are both 
performed. In any case, the large uncertainties on and make it hard to determine the 
magnitude of the redshift evolution of Pth,c- 

The thermal pressure can be used to predict the amount of cosmic microwave 
background (CMB) diminution, also known as the Sunyaev-Zel'dovich (S-Z) effect, that is 
expected from the cluster. The reduction in the cosmic radiation in the direction of 

a cluster is given by {e.g. Sunyaev & Zel'dovich 1980; Sarazin 1988; Rephaeli 1995) 

A7,/7, = G{u;)y. (13) 

Here Ii, is the specific intensity of the CMB at the observing frequency y is the 
Comptonization parameter, and the function G{uj) is 
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e'^ - 1 

where u = hu/kTr, and Tr = 2.73 K is the CMB temperature (Mather et al. 1990). 
The Compton y parameter is given by 

y = / -^^arUedl = ( — — — ) / Pthdl. (15) 

J meC^ rUeC^ \ne + Ui J J 

where gt — (87r/3)[e^/(meC^)]^ is the Thompson electron scattering cross section, T is the 

cluster gas temperature, is the electron density, and rii is the ion density. Using the 

thermal pressure in a cluster with an isothermal- King gas distribution (eq. [11]), this gives 

/ Pthc \ f rc \ r(3/3o/2 - 1/2) 

t/«6.3xlO-M — ]-—^ ^(i + ^2)-3/.o/2+i/2^ ^;L6a) 

V 10-1° dyne cm-2 / Vo.25Mpc/ T{3(3o/2) 

where x = r/rc, and the gas is taken to have solar abundance so that Ue = 1.21nj. For a 

typical value of /3o = 2/3, 

-4 / Pth,c \ I ^r. 



y « 1.12 X 10-^ . , (1 + x2)-V2. (16b) 

^ \lO-iodynecm-V V0.25Mpc^' ^ ^ ^ ^ 
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At low frequency, clusters with strong radio sources are generally avoided for 
measurements of the Sunyaev-Zel'dovich effect because emission from the radio source 
masks the microwave diminution. Such contaminations from radio sources are reduced at 
high frequency since emission from steep-spectrum radio sources, such as the FRII sources 
studied here, decreases rapidly with increasing frequency The Sunyaev-Zel'dovich Infrared 
Experiment (SuZIE) can measure the S-Z effect around 140 GHz, where the amount of CMB 
intensity diminution is near its peak. At this frequency, a cluster with an isothermal- King 
density profile will cause a CMB intensity diminution of about 

Pth,c 



AI^ 7.0 X 10"^^erg s-^Hz-^cm-^sr"^ 

r(3/?o/2 - 1/2) 2^-3/3o/2+l/2 

r(3/?o/2) ^'^""^ 
and for a typical value of /3o = 2/3, 



10-^° dyne cm-2 J \0.25Mpc 



;i7) 



A7, ^ 1.24 X 10-^«erg s-^Hz-^cm-^sr-i ( -fi^ ( (1 + x^)-^^. 

^ V 10-10 dyne cm" V V 0-25 Mpcy ^ ^ 

(18) 

For a detector with a FWHM beam size 9b, the total amount of CMB diminution 
within the beam, defined as AFi,, is roughly 

9b \^ f Pth,c 



AF,f«/-24mJy ^ , (19) 

Vl-7 arcminy VIQ-^" dyne cm" V V 0-25 Mpc/ ' ^ ^ 

where / is the beam dilution factor, defined as the ratio of the average AI^, within the 
beam to the peak value at the cluster center; and the SuZIE FWHM beam size at 140 GHz 
of 1.7 arcmin {e.g. Holzapfel et al. 1997) is used to calculate the numerical value. For a 
Gaussian beam, when the HWHM beam size corresponds to Ir^ f ~ 0.7, and when the 
HWHM beam size corresponds to 2rc, / ^ 0.5 (see Rephaeh 1987). At z < 2, the SuZIE 
beam radius is less than ~ 320/i-ikpc for qo = with no cosmological constant. Thus a 
cluster with Pth,c'^ 10~i°dyne cm~^ and 0.25 Mpc will have a AF^ of about (15 to 
20) mJy within the SuZIE beam. The clusters surrounding the FRII sources in our sample 



-23- 



have Pth,c on the order of 10~^° /i^/^ dyne cm~^ and the data also suggest that Tc increases 
with redshift, reaching about 300/i~^kpc at 2;~2. The expected CMB diminution within 
the SuZIE beam at 140 GHz for these clusters ranges from several to tens of mJy. These 
clusters make good candidates for SuZIE observations if the fluxes from the radio sources 
at 140 GHz or their uncertainties are small compared with the expected S-Z effect signals. 

We are currently in the process of searching for high frequency data on the radio 
sources in our sample in oder to identify possible SuZIE observation candidates. One likely 
candidate that comes from a preliminary search in the published literature is 3C239. The 
expected CMB diminution for the cluster surrounding it is about 30 mJy within the 
SuZIE beam at 140 GHz. In contrast, extrapolation of the observed spectrum of 3C239 as 
given by observations at 178 MHz, 10.7 GHz, and 14.9 GHz (Kellermann & Pauhny-Toth 
1973; Genzel et al. 1976; Laing, Riley, and Longair 1983), gives a 140 GHz flux of only 6 
mJy. This is likely to be an overestimate since synchrotron aging and inverse Compton 
cooling can both cause the high frequency spectral index to become steeper than that at 
low frequency. The effects of inverse Compton cooling can be especially important since this 
source is at high redshift {z — 1.79), where the energy density of the microwave background 
is much higher than at low redshift. Thus it appears that emission from this radio source 
is weak compared with the expected CMB diminution. Observations of the radio source at 
more frequencies above 14.9 GHz can help to better constrain its 140 GHz flux. 

5. Gravitational Mass of the Host Cluster 
5.1. Theory 

The powerful classical double radio sources studied here are in cluster-like gaseous 
environments (see WDW97a, b). The studies of these sources provide information on 
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density and temperature of the ambient gas. Thus it is possible to estimate one of the most 
important parameters of the cluster, the total gravitational mass, including dark matter. 

The total mass can be estimated using the density and temperature profile of the 
intracluster medium (ICM) if the gas is in hydrostatic equilibrium. The sound crossing 
time in the ICM is usually short compared to the age of a high-redshift cluster, and the 
morphology of the X-ray emission from many low-redshift clusters is often smooth. Thus 
it is generally believed that hydrostatic equilibrium is a good approximation of the state of 
the ICM for many clusters that are not cooling flow clusters. Assuming spherical symmetry, 
hydrostatic equilibrium requires 

1 dP _ _GMt{r) 
Pg dr ' 

where Pg is the gas density and Mt(r) is the total gravitating mass within r. This means 
that 

J^T(n^\ rd\ogpg{r) ^ d\ogT{ry ^^^^ 



M .X _ kT{r) 

* GpiTLp 



d log r d log r 

where k is the Boltzman constant, p is the mean molecular weight of the gas in amu, rUp is 
the proton mass, and T(r) and pg{r) are the temperature and density profiles of the cluster 
gas. 

The most commonly used hydrostatic model of the ICM is the isothermal (3- 
model (Cavaliere & Fusco-Femiano, 1976,1978; Sarazin & Bahcall 1977), where the 
clusters gas is isothermal and the density profile follows a modified King model, i.e.. 
Ha — nc\\ + (r/rc)^]~^/^^°. Here ric is the core density and Tc is core radius of the gas 
distribution. Using this model for the ICM, eq. (21) becomes 

Results from numerical simulations of cluster formation suggest that departures of the ICM 
from hydrostatic equilibrium are usually small, and the mass estimated using the standard 



-25- 



/3-model is rather accurate {e.g. Navarro, Prenk & White 1995; Schindler 1996; Evrard, 
Metzler, & Navarro 1996). Note however, significant temperature dechne is observed to 
occur at outer regions of clusters {e.g. Markevitch et al. 1998). Thus, we only use eq. (22) 
to estimate the cluster mass within r rather than extrapolating to large radii. 

The mass estimated using eq. (22) is not accurate if the cluster is cooling at the point 
where the cluster temperature is measured. Inside a cooling fiow region, the temperature 
measured does not indicate the gravitational potential, since hydrostatic equilibrium 
conditions do not apply. Considering cooling only by thermal bremsstrahlung, which 
dominates other mechanisms at typical cluster temperatures, the cooling time of the cluster 
gas can be estimated using 

*w-2.1xlO-(^)"'(-^)-'Gy, (23) 

where T and Ua are the ambient gas temperature and electron density, respectively (see 
WDW97a). Note this equation does not depend on the Bubble's constant, provided that 
the estimates of and T do not depend on Hubble's constant. 

Given the ambient gas temperature and density estimates, the cooling time for the 
sources in our sample can be estimated (see WDW97a). A cooling time less than the age 
of the universe at that redshift means that the FRII source is in a cooling flow region. To 
estimate the age of the universe at a given redshift, an open empty universe is assumed, 
with a current age of 14 Gyr, which corresponds to /i ~ 0.7 for = with zero cosmological 
constant. 

Most sources in our sample do not appear to be cooling at the position where the 
temperature is measured. Thus their mass estimates are likely to be reliable. A few sources, 
including the low-redshift source Cygnus A, appear to lie in cooling flow regions. 

Figure 9 shows the total cluster mass, including dark matter, within radius r, Mt(r), as 
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a function of the core-hot spot separation r. Since the sources appear to he at the centers 
of clusters of galaxies (Daly 1995; WDW97a,b), the core-hot spot separation is used to 
estimate the distance from the cluster center. The radio spectral index is redshift- corrected, 
and the cluster core radius is taken to be Tc — 50(1 -|- zY-^ kpc (see §4.2). Clusters that 
are cooling at the position where the temperature is measured are marked with pentagons. 
The plots of Mt{r) vs. r for other models, either with or without an a—z correction and/or 
different core radius evolution, are presented in Figures 10 through 12. 

The cluster mass increases with r for all models. The increase is more rapid for models 
with an evolving core radius than models with a fixed core radius. This merely reflects the 
fact that with an increasing Tc, more sources lie within the core region, where the increase 
of mass with radius is rapid, and varies roughly as M oc r^. 

To determine the redshift evolution of the cluster mass, a three parameter fit of 
Mt{r) oc r"i(l + zY^'^P'^f^ is performed since the cluster mass within r, Mt(r), is a function 
of r, and may also be affected by radio power selection effects. Results are listed in Table 2 
for all models considered. A two parameter fit of Mt{r) oc r"i(l -|- zY'^ is also performed, 
and results are hsted in Table 3. For completeness, fits including all sources, all sources 
except Cygnus A, and all sources that are not in cooling flow regions, are all listed in 
the table. Note that Cygnus A appears to be in a cooling flow region for all the models 
considered. Results obtained excluding cooling flow clusters are probably the correct ones 
to consider. 

The sample of clusters with mass estimates is rather small, especially when only 
non-coohng flow clusters are considered. Thus it is not surprising to see that the best-flt 
values of n2 have rather large uncertainties, which makes it hard to draw definitive 
conclusions about the redshift evolution of the cluster mass. Note though, the current 
data do not indicate any negative evolution of the cluster mass out to a redshift of about 
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two. This is consistent with results obtained by this group (Daly 1994; Guerra & Daly 
1996, 1998; Guerra, Daly, & Wan 2000; and Daly, Guerra, & Wan 1998), who study the 
characteristic size of powerful classical double radio sources and find the data suggest a low 
value of ftjn] a universe with = 1 is ruled out at 99 % confidence (see Guerra, Daly, & 
Wan 2000). 

There are some indications from the data that the redshift evolution of the cluster core 
radius, and the redshift-correction on the radio spectral index are favored. The gas mass 
within r, defined as Mg{r), can be estimated using 

Mg{r) = inpcrlix - tan-^(x)], (24) 

where Pc is the central gas density, Tc is the core radius, x = r/rc, and a typical value of 
/3o = 2/3 is used for the King density profile. Knowing the gas mass and the total mass 
within r, the gas mass fraction within r can be estimated. The gas mass fraction within 
radius r as a function of redshift is shown in Figure 13, where the radio spectral index is 
redshift-corrected and the core radius is taken to be Vc — 50(1 + z)^-^ kpc. 

The values of gas fraction shown in the figure are consistent with observed values for 
inner regions of many clusters {e.g. Donahue 1996), whereas those for models without a 
redshift-correction on the radio spectral index, and without a core radius evolution seem to 
be low compared with observed values. 

Several factors may contribute to the slow decrease of gas fraction with redshift that 
is seen in Figure 13. First, the increase of cluster core radius with redshift means that 
the gas fraction estimated for the high-redshift clusters is over a larger fraction of the 
cluster core than for the low-redshift clusters. The gas fraction estimated at high redshift 
mainly represents the gas fraction inside the cluster core, and it is known that gas fraction 
increases with increasing radius in clusters {e.g. David, Jones, & Forman 1995). Second, 
the mass contribution from galaxies is higher at the cluster center than at large radii {e.g. 
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Loewenstein & Mushotzky 1996). Thus by sampling a larger fraction of the cluster core at 
higher redshift, a larger fraction of the total baryon mass is not taken into account at higher 
redshift. Further, the decrease in gas fraction with redshift could be due to the stripping 
of gas from galaxies in the cluster core over time. Finally, the cluster gas becomes more 
concentrated toward the cluster center as it cools, which also causes the gas fraction in the 
cluster center to increase with time. 

These results on cluster mass and gas fraction obtained above are still preliminary 
due to the small size of the sample. More sources with estimates on the ambient gas 
temperature, and hence cluster mass, will help to test these results. 

6. Summciry and Discussion 

Several key parameters of an FRII source and its gaseous environment are studied in 
this paper. 

Direct estimates of the beam power, which measures the energy extraction rate from 
the AGN by the jet, are obtained. Typical beam powers of about 10^^ h""^ erg are found 
for the sources in the sample. No strong correlation is seen between the beam power and the 
linear size of the source, which is consistent with the beam power being roughly constant 
throughout the lifetime of a source. 

There is a trend for the beam power to increase with redshift, which is significant even 
after excluding radio power selection effects. The magnitude of this redshift evolution, 
however, is not well determined by the current data. It is well known that the quasar 
luminosity function undergoes strong evolution between z ^ and z ^ 2, with the 
high-redshift quasars being more luminous than their low-redshift counterparts {e.g. 
Schmidt & Green 1983; Yee & EUingson 1993; La Franca & Cristiani 1997). This suggests 



-29- 



that high-redshift AGNs are more powerful than low redshift ones. Thus it is perhaps not 
surprising to find that the beam power, or energy per unit time channeled into the jet by 
the AGN, is also higher at high redshift. 

The relationship between the beam power and the radio power is not well constrained 
after their correlations with redshift are taken into account. The two parameter fit of 
Lj (X (1 + -z)"^-Pi"8 yields values of Up with large uncertainties. Thus it is not clear at 
present how the beam power and the radio power are related. However, the large amount 
of scatter seen in the beam power-radio power relation suggests that radio power is not an 
accurate measure of beam power. 

The beam power Lj can be used to estimate the total lifetime of an outflow produced 
by an AGN. Following Daly (1994) and Guerra & Daly (1996, 1998), the total time for 
which the outflow will occur, i*, is related to Lj: oc Lj'^*^^, where is estimated to be 
about 1.75 ± 0.25 (Guerra, Daly, & Wan 2000). Typical lifetimes of about 10^ to 10^ years 
are obtained for the sources studied here. This would is almost precisely the same lifetime 
as would be obtained by dividing the average size of all FRIIb sources at similar redshift by 
the rate of growth of the source under consideration. The hfetime of the outflow decreases 
with redshift, which explains the decrease of the average size of powerful extended SCR 
sources with redshift. The sources at high redshift produces more powerful jets for a shorter 
period of time, which results in smaller average sizes than low-redshift sources. 

A new method of estimating the thermal pressure of the ambient gas in the vicinity of 
powerful classical double radio source using only single frequency radio data is presented. 
The pressure is given by the product of the ambient gas density, estimated using ram 
pressure conflnement of the radio lobe, and the ambient gas temperature, estimated using 
the Mach number for the source and the lobe advance speed. It turns out that the 
lobe propagation velocity cancels out of the product UaT, and the ambient gas 
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pressure can be estimated by studying the shape of the radio bridge, and the 
non-thermal pressure of the radio lobe (see §4). Thus, the thermal pressure, the 
product of the ambient gas density and temperature, depends only on the non-thermal 
pressure in the radio lobe and the geometrically determined Mach number of lobe advance, 
both of which can be estimated using single frequency radio data. This new method to 
estimate the thermal pressure does not require a spectral aging analysis, and provides a 
powerful tool to probe the environments of FRII sources. The thermal pressure estimated 
for Cygnus A using this new method agrees with that obtained using X-ray data (see §4.2.). 

Thermal pressures on the order of 10~^° /i^/^dyne cm~^, typical of gas in low-redshift 
clusters of galaxies, are found for the gaseous environments of the FRII sources studied 
here. There are hints from the current data that the gradient of the composite pressure 
profile is less steep at high redshift than at low redshift, which can be explained by an 
increase of the cluster core radius with redshift. The current data are consistent with a 
core radius evolution from ~50/i~^kpc at 2;~0 to ~250/i~^kpc at 2;~2, which agrees with 
results obtained from studies of the ambient gas density (WDW97b). WDW97b find that 
the data can be described by a model where the core gas density decreases and the core 
radius increases so that the core gas mass remains roughly constant. 

The thermal pressures obtained here can be used to estimate the amount of CMB 
diminution expected from the clusters surrounding the FRII sources in the sample. 
Contaminations from the radio sources can be reduced by observing at high frequency, 
such as 140 GHz, since emission from the radio sources decreases rapidly with increasing 
frequency. The cluster surrounding a source in our sample can be detected by SuZIE 
observations at 140 GHz if the flux from the radio source at this frequency or its uncertainty 
is small compared with the expected S-Z effect signal. A search for high-frequency data 
of the sources in our sample is ongoing in oder to identify possible SuZIE observation 
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candidates. Preliminary results suggest that 3C239 is a good candidate. 

The gravitational or total mass of the surrounding cluster is estimated for the sources 
in the sample, assuming hydrostatic equilibrium conditions for the gas. Masses of up to 
IO^^Mq are found for the central regions (r < 250 h^^ kpc) of the clusters, consistent with 
typical values for low-redshift clusters. The redshift evolution of the cluster mass is not 
well determined. Current data do not indicate negative evolution of the cluster mass. This 
is consistent with results obtained by members of this group (Guerra & Daly 1996, 1998; 
Guerra, Daly, & Wan 2000; Daly, Guerra, & Wan 1998), who study the characteristic size 
of FRII sources and find the data strongly favor a low value for the density parameter fl^; 
a universe with = 1 is ruled out at 99 % confidence. Note, however, that with the study 
presented here, we only study the central regions of clusters. 

The values of gas mass fraction obtained for the clusters surrounding the FRII sources 
studied here are consistent with observed values for central regions of clusters, after the 
correlation between the radio spectral index and redshift is taken into account, and a 
redshift evolution of the cluster core radius is considered. This suggests that the cluster 
core radius evolution and the effects of the radio spectral index-redshift correlation are 
important. 

The gas mass fraction seems to decrease slowly with redshift. The increase of cluster 
core radius with redshift can be one cause. Another possible explanation is that a large 
fraction of the gas is still in the galaxies at high redshift, and is later stripped from the 
galaxies into the cluster at low redshift. Coohng of the cluster gas also tends to cause the 
gas mass fraction in the cluster center to increase with time. 

The results on mass and gas fraction should be taken as preliminary because of the 
small size of the sample. 
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Table 1 — Continued 



Source 
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^OgPth" 


LogPth,c(A)d 














45.15 ±0.15 


45.02 ±0.14 


-11.18 ±0.28 


-10.30 ±0.28 


-10.76 ± 0.28 


3.0 ±0.6 


3.5 ±0.7 


3C 334 


Q 


0.555 


44.94 ±0.12 


44.81 ±0.12 








4.0 ±0.5 


4.7 ±0.6 








45.05 ± 0.22 


44.89 ± 0.20 








3.4 ±0.9 


4.2 ± 1.1 


3C 356 


G 


1.079 


45.62 ±0.21 


45.41 ±0.21 


-10.51 ±0.22 


-9.14 ±0.22 


-10.01 ±0.22 


1.6 ±0.4 


2.1 ±0.6 








45.88 ±0.21 


45.69 ±0.21 


-10.59 ±0.17 


-9.54 ±0.17 


-10.29 ±0.17 


1.1 ±0.3 


1.4 ±0.4 


3C 405 


G 


0.056 


45.07 ±0.13 


45.06 ±0.13 


-10.31 ±0.19 


-10.07 ±0.19 


-10.10 ±0.19 


3.3 ±0.6 


3.4 ±0.6 








44.99 ±0.12 


44.97 ±0.12 


-10.04 ±0.17 


-9.74 ±0.17 


-9.78 ±0.17 


3.7 ±0.6 


3.8 ±0.6 


3C 427.1 


G 


0.572 


44.72 ±0.13 


44.60 ±0.13 


< -9.81 


< -9.52 


< -9.72 


5.3 ±0.9 


6.2 ± 1.0 








44.83 ±0.11 


44.71 ±0.11 


-10.37 ±0.09 


-10.05 ±0.09 


-10.27 ±0.09 


4.6 ±0.6 


5.4 ±0.7 



^Logarithm of the luminosity in directed kinetic energy of the jet in h ^ erg s ^. 

''Logarithm of the luminosity in directed kinetic energy of the jet in h^^ erg s~', with redshift-corrected radio spectral indices 
(see §3.1). 

''Logarithm of the thermal pressure of the ambient gas in the vicinity of the radio lobe in h^/"^ dyne cm~^. The 3a upper 
bounds are listed for sources without detections. 

Logarithm of the thermal pressure at the center of the surrounding cluster in h^/'^ dyne cm~^, assuming a cluster core radius 

of 50/1-1 ]^p^ (ggp p 2). 

''Logarithm of the thermal pressure at the center of the surrounding cluster in /i'*/'^ dyne cm" ^, assuming a cluster core radius 
of 50 (1 + 2)1-6 /i-i kpc (see §4.2). 

^Characteristic lifetime of the outflow in (1 ± 0.25) x lO^h^/^ yrs, where the uncertainty is that on the normalization factor 

C. 

^Characteristic lifetime of the outflow in (1 ± 0.25) x lO'^h'^/^ yrs, with redshift-corrected radio spectral indices (see §3.1). 
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Table 2. Results Obtained by Fitting Mt{r) oc r"i(l + z)''^P\ 



type"* No> a-z?" Tc-zI'^ nx" n-i^ n^^ Xr 



N.C.Cl. 11 Yes Yes 2.38 ± 0.70(0.52) 6.86 ± 5.48(4.1) -1.31 ± 1.57(1.17) 0.56 
G+Q-C 14 Yes Yes 2.95 ± 0.31(0.23) 2.72 ± 3.80(2.84) -0.08 ± 1.01(0.76) 0.56 
G-l-Q 16 Yes Yes 2.88 ± 0.30(0.23) -0.73 ± 0.96(0.73) 0.80 ± 0.37(0.28) 0.58 



N.C.Cl. 13 No Yes 3.05 ± 0.36(0.28) 
G+Q-C 14 No Yes 2.95 ± 0.31(0.23) 
G+Q 16 No Yes 2.89 ± 0.30(0.23) 



4.46 ± 4.11(3.18) -0.13 ± 1.06(0.82) 0.60 
3.66 ± 3.82(2.88) 0.02 ± 1.02(0.77) 0.57 
0.49 ± 0.96(0.73) 0.83 ± 0.37(0.28) 0.58 



N.C.Cl. 


11 


Yes 


No 


1.36 ± 0.63(0.52) 


7.40 ± 4.55(3.72) 


-0.93 ± 1.33(1.09) 


0.67 


G+Q-C 


14 


Yes 


No 


2.12 ± 0.27(0.23) 


3.18 ± 3.33(2.85) 


0.44 ± 0.90(0.77) 


0.73 


G+Q 
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Yes 


No 


2.11 ± 0.26(0.21) 


1.63 ± 0.91(0.75) 


0.84 ± 0.33(0.27) 


0.68 



N.C.Cl. 


13 


No 


No 


2.17 ± 0.31(0.28) 


4.45 ± 3.40(3.10) 


0.47 ± 0.88(0.80) 


0.83 


G+Q-C 


14 


No 


No 


2.12 ± 0.27(0.23) 


4.09 ± 3.20(2.77) 


0.53 ± 0.86(0.74) 


0.75 


G+Q 


16 


No 


No 


2.11 ± 0.27(0.22) 


2.84 ± 0.91(0.76) 


0.86 ± 0.33(0.27) 


0.69 



''Type of sources included in the fit: "G+Q" refers to all sources with temperature estimates, including 
radio galaxies and radio-loud quasars, "G+Q-C" refers to sources other than Cygnus A, and "N.C.Cl." refers 
to clusters that are not cooling at the position where the temperature is measured. 

''Number of data points used for the fit. 

'^Whether the rcdshift-corrcction on the radio spectral index is applied. 

■^Whether a rcdshift evolution of the cluster core radius Tc is considered. The core radius is taken to be 
50 (1 + z)^'^h^^ kpc if an evolution is considered, and 50 h^^ kpc if no evolution is considered (see §4.2). 

and its error. The number in parenthesis is the error on rii times -^/xF which includes the effect of the 
reduced x^j defined as Xr: of the fit. 

^'^same as note c for n2 and ns, respectively. 
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Table 3. Results Obtained by Fitting Mt{r) oc (1 + 



type^ No.*^ a-z'?'^ Vc-zl'^ ni^ Xr 



N.C.Cl. 11 Yes Yes 2.70 ± 0.59(0.45) 2.39 ± 1.24(0.94) 0.58 
G+Q-C 14 Yes Yes 2.95 ± 0.31(0.22) 2.43 ± 1.05(0.75) 0.51 
G+Q 16 Yes Yes 2.67 ± 0.29(0.28) 0.63 ± 0.73(0.69) 0.90 



N.C.Cl. 13 No 
G+Q-C 14 No 
G+Q 16 No 



Yes 3.04 ± 0.35(0.26) 
Yes 2.95 ± 0.31(0.22) 
Yes 2.67 ± 0.29(0.28) 



3.99 ± 1.16(0.86) 0.55 
3.74 ± 1.06(0.76) 0.52 
1.90 ± 0.73(0.70) 0.93 



N.C.Cl. 11 Yes 
G+Q-C 14 Yes 

G+Q 16 Yes 



No 1.60 ± 0.52(0.42) 
No 2.13 ± 0.27(0.22) 

No 1.88 ± 0.25(0.27) 



4.30 ± 1.04(0.84) 0.65 
4.75 ± 0.83(0.68) 0.68 
3.34 ± 0.62(0.66) 1.13 



N.C.Cl. 


13 


No 


No 


2.20 ± 0.31(0.27) 


6.19 ± 0.94(0.82) 


0.77 


G+Q-C 
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No 


No 


2.13 ± 0.27(0.23) 


6.00 ± 0.83(0.70) 


0.72 


G+Q 


16 


No 


No 


1.88 ± 0.25(0.27) 


4.57 ± 0.62(0.67) 


1.16 



^ -^same as those in Table 2. 
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Fig. 1. — The luminosity in directed kinetic energy of the jet vs. the core-hot spot separation. 
The best-fitted lines for all galaxies, all galaxies except Cygnus A, and all radio-loud quasars, 
are shown with a sohd line, a dashed line and a dotted line, respectively. Slopes of the best- 
fitted lines are labeled. The number in parenthesis is the uncertainty on the slope times 
-^xf, which includes the effect of a reduced x^, defined as x^, greater than one for the fit. 
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Fig. 3. — The luminosity in directed kinetic energy of the jet vs. the radio power at 178 
MHz. 
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Fig. 4. — Total lifetime of the coUimated outflows vs. redshift. 
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Fig. 5. — Thermal pressure of the ambient gas in the vicinity of the radio lobe vs. the core- 
hot spot separation. The best-fitted lines for all sources, including galaxies and quasars, and 
all sources except Cygnus A are shown with a dashed line and a dash-dot line, respectively. 
Only sources with detections of the Mach number M, and hence Pth, are included in the fit. 
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Fig. 6. — Thermal pressure of the ambient gas in the vicinity of the radio lobe vs. the 
core-hot spot separation for two redshift bins {z < 0.9, z > 0.9). In the low-redshift bin, 
slopes of the best-fitted lines are —1.43 ±0.21 (0.27) with a reduced of 1.58 for all sources, 
and —1.35 ± 0.25(0.32) with a reduced of 1-65 for all sources except Cygnus A. In the 
high redshift bin, the slope of the best-fitted line is —0.41 ± 0.20(0.30) with a reduced of 
2.29. Only sources with detections of Pth are included in the fit. 
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Fig. 7. — Thermal pressure of the ambient gas at the cluster center vs. redshift for an 
isothermal gas distribution that follows the King density profile. A non-evolving cluster 
radius of bOh"^ kpc, and Po — 2/3 are used. The same plots for an evolving core radius 
of 50(1 + zY'^h'^ kpc are very similar to these and are not shown here. Bounds are not 
included in the fits. 
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Fig. 8. — Thermal pressure of the ambient gas at the cluster center vs. the radio power at 
178 MHz for an isothermal gas distribution that follows the King density profile. A non- 
evolving cluster radius of 50h~^ kpc, and Po — 2/3 are used. The same plots for an evolving 
core radius of 50(1 + zy^^h~^ kpc are very similar to these and are not shown here. Bounds 
are not included in the fits. 
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Fig. 9. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas 
distribution that follows the King density profile with an evolving cluster core radius of 
50 (1 + kpc and /Sq = 2/3 is used. The radio spectral index is redshift-corrected. 

Sources that are cooling at the point where the cluster temperature T is measured are 
marked with pentagons. The best-fitted lines for all sources, including both galaxies and 
quasars, all sources except Cygnus A, and all sources that are not cooling, are shown with 
a dashed line, a dot-dash line, and a solid line, respectively. Only sources with detections of 
T, hence mass, are included in the fit. 
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Fig. 10. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas 
distribution that follows the King density profile with an evolving cluster core radius of 
50 (1 -I- zy-^h~^ kpc and Pq — 2/3 is used. The radio spectral index is not redshift-corrected. 
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Fig. 11. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas 
distribution that follows the King density profile with a fixed cluster core radius of 50 kpc 
and /So = 2/3 is used. The radio spectral index is redshift-corrected. 
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Fig. 12. — Cluster mass within r vs. the core-hot spot separation r. An isothermal gas 
distribution that follows the King density profile with a fixed cluster core radius of 50 kpc 
and /So = 2/3 is used. The radio spectral index is not redshift-corrected. 
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Fig. 13. — Gas fraction {Mgas/Mt) within radius r vs. redshift. The radio spectral index is 
redshift-corrected, and an isothermal gas distribution that follows the King density profile 
with an evolving cluster core radius of 50 (1 + zY'^h^^ kpc and Pq = 2/3 is used. 



